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I.  INTRODUCTION 

\h  e  theory  of  optimum  arrays  became  widely  known  in  the  underwater  acoustics 
community  through  the  work  of  Bryn  and  Mermoz  [i7].  The  relationship  of  their  work 
to  other  detection  and  estimation  problems  has  been  discussed  in  l/f\.  Optimum  array 
processing  structures  use  detailed  information  about  the  signal  and  noise  fields.  Since 
this  information  is  not  known  precisely  in  advance,  one  is  led  naturally  to  adaptive  beam- 
formers  which  continually  adjust  their  parameters  based  upon  on-line  measurements  of 
some  kind.  Since  adaptive  processors  are  continually  adjusting,  it  is  natural  to  question 
how  sensitive  performance  is  to  small  variations  of  die  signal  field,  noise  field  and  sys¬ 
tem  parameters  from  their  assumed  or  estimated  values. 

The  question  of  sensitivities  has  been  examined  in  the  past  in 

conjunction  with  *  super  -directive*  arrays.  An  attempt  will  be  made  to  point  out  the  re¬ 
lationship  of  the  results  of  this  paper  to  those  earlier  results.  The  emphasis  in  this  paper 
is  on  receiving  arrays. 

In  Section  II  an  introduction  to  the  problem  is  provided  using  an  intuitive  approach. 
The  performance  measures  of  array  gain  and  output  power  are  used. 

Section  III  presents  the  principal  sensitivity  results.  The  approach  is  to  take  par¬ 
tial  derivatives  of  the  gain  and  output  signal  power  with  respect  to  the  size  of  signal,  noise 
and  steering  perturbations. 

Section  IV  discusses  the  problem  of  signal  suppression  which  arises  in  passive 
adaptive  processors  when  measurements  of  signal-plus -noise  are  used  when  noise  only 
measurements  are  desired.  Interference  rejection  is  also  discussed. 

A  number  of  optimization  problems  are  discussed  in  Section  V.  Particular  em¬ 
phasis  is  given  to  formulations  which  somehow  take  sensitivity  into  aqfco\int.  It  is  pointed 
out  that  the  simple  beamformer  structure  is  no  longer  optimum  in  most  o!^  these 
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problems.  A  more  general  array  structure  is  suggested  which  possesses  optimality 
properties. 

Vector-matrix  notation  similar  to  that  in  [7]  will  be  used  throughout  this  paper. 
Column  vectors,  i.  e. ,  matrices  with  only  one  column  are  designated  by  underlined  lower 
case  letters  such  as  k  and  m.  Other  matrices  are  designated  by  underlined  upper  case 
letters  such  as  P,  Q  and  R.  The  asterisk  is  used  to  denote  complex  conjugate  transposi¬ 
tion.  Thus,  (k*k)  is  a  scalar  and  (mm*) is  a  dyadic  square  matrix. 

n.  PRELIMINARY  DISCUSSION 

The  receiving  array  consists  of  an  arbitrary  arrangement  of  M  sensors.  When  a 
signal  is  present  the  waveform  at  the  output  of  the  i-th  sensor  is 
*j(t)  =  v.(t)  +  n.(t) 

where  v.(t)  is  the  signal  component  and  n.(t)  is  the  noise  component.  The  noise  component 
includes  both  the  effects  of  the  external  noise  field  and  internal  electronic  noise.  The  out¬ 
puts  of  all  M  sensors  may  be  represented  by  the  following  vector  equation: 

x(t)  =  v(t)+ n(t)  (1) 

The  detection  problem  is  usually  formulated  as  deciding  whether  or  not  the  signal  com¬ 
ponent  v(t)  is  present.  Estimation  problems  arise  when  one  seeks  to  obtain  information 
about  particular  parameters  of  v(t).  The  relationships  among  a  number  of  detection  and 
estimation  problems  are  discussed  in  [7].  There  it  is  shown  that  the  same  beamforming 
processor  is  optimum  for  a  variety  of  different  detection  and  estimation  problems. 

The  signal  vector  v^(t)  and  the  noise  vector  n(t)  are  assumed  to  have  mean  values  of 
zero  and  to  be  statistically  independent  of  each  other.  Because  the  analysis  of  array  pro¬ 
cessors  is  sometimes  simplified  by  working  in  die  frequency  domain,  it  is  convenient  to 
work  with  the  cross -spectral  density  matrices  of  v(t),  n(t)  and  x(t)  which  are  denoted 
P(ui),  Q(u)  and  R  («j)  respectively.  Then 

R(u)  =  P(u>)  +  Q  (uj)  (2) 

when  the  signal  is  present.  P(u)  and  Q(u)  may  be  written  as 

P(u)  =  cr*(u>)  P(u>)  (3) 

"  8  “ 

Q(w)  =  o^(u>)  Q(w)  (4) 

2 

where  a  (u)  is  the  noise  power  spectral  density  averaged  over  die  M  sensors  so  that 
n  2  ~ 

a  (u)  =  trace  Q(u)/M  (5) 

2  n 

and  <J  (w)  is  similarly  defined  as  the  signal  power  (or  energy)  spectral  density  averaged 

S 
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over  the  M  sensors.  Thus,  P(u>)  and  Q(u>)  are  normalized  to  have  their  traces  equal  to  M, 

2  2 

the  number  of  sensors  in  the  array.  The  quantity  a  (u)/a  (u>)  is  the  input  signal -to- 

s  n 

noise  spectral  ratio  (S/N).  . 

A  device  which  has  the  configuration  shown  in  Figure  1  will  be  called  a  beamfbrmer. 
In  general,  it  consists  of  a  set  of  filters,  one  for  each  sensor,  followed  by  a  summation 
device.  The  transfer  function  of  the  filter  on  the  i-th  sensor  is  k*  (u>).  Two  quantities  of 
prime  interest  are  the  output  spectral  density 

z(u)  =k*(u>)  R(u)  k  (u)  (6) 

and  the  array  gain 


G(u) = • 


k*(u>)  P(w)  k(u) 
k*(w)  Q(u)k(u) 


where,  from  Figure  1,  the  steering  vector  k*(u)  is  the  row  vector  of  filter  transfer  func¬ 
tions  jk*(u>),  *•••»  kj^(u>)  |  .  The  sensitivity  analysis  will  be  concerned  with  determining 
the  effects  on  z(u)  and  G(u)  of  varying  P(j),  Q(u)  and  k(u)  from  their  nominal  values. 

The  array  gain  given  by  (7)  is  the  ratio  of  the  output  signal-to -noise  spectral  ratio 
to  the  input  signal-to-noise  spectral  ratio.  It  is  easily  seen  to  be  invariant  if  k*(u)  is 
multiplied  by  a  scale  factor.  In  considering  sensitivity  problems  it  is  sometimes  useful  to 
normalize  both  the  numerator  and  denominator  of  (7)  by  dividing  by  the  magnitude  squared 
of  k*(u).  Then  (7)  becomes 

k*(u)  P(u>)  k(u) / (k*(u)  k(u) ) 

°<W)  =  k*(u)Q(u)k(u) /(k*(u;)  k(u)) )  (8) 

The  numerator  of  (8)  is  the  normalized  signal  response  and  the  denominator  is  the  normal¬ 
ized  noise  response.  We  shall  see  that  the  numerator  and  denominator  of  (8)  individually 
or  more  specifically  their  reciprocals  are  of  fundamental  importance  in  sensitivity  analy¬ 
sis. 

One  of  the  reasons  for  working  in  the  frequency  domain  is  that  the  signal  vector  v(t) 
is  usually  assumed  to  be  related  to  a  scalar  signal  s  (t)  at  some  source  point  in  space  by  an 
equation  of  the  following  form: 

v(t)  =  y*m(t  -  t)  s(t)  dr  (9) 

where  m^t)  accounts  for  the  propagation  from  the  source  to  the  i-th  sensor  and  the  re¬ 
sponse  of  the  i-th  sensor  itself.  In  the  ideal  case  of  non-dispersive  propagation  and  dis¬ 
tortion-free  sensors,  m^t)  is  a  simple  time  delay  6  (t  -  T^).  Whenever  v(t)  is  related  to 
a  scalar  s(t)  by  a  known  transformation  m(t)  as  in  (9)  P(uj)  is  simply  the  dyadic 
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P  (u)  =  m(u>)  m*(u))  (10) 

where  m(u)  is  the  Fourier  transform  of  m(t)  normalized  so  that  m*(u)m(u))  =  M. 

The  noise  field  may  be  the  sum  of  a  number  of  components.  Some  components  of 
particular  significance  are  spatially  white  noise,  isotropic  noise,  and  point  source  noise. 

Spatially  white  noise  is  uncorrelated  from  sensor  to  sensor  and  of  equal  intensity  at 
each  sensor.  It  has  the  identity  matrix  for  its  normalized  cross -spectral  density  matrix. 
The  subscript  w  on  =_I  will  be  used  to  denote  this  type  of  noise  component. 

Another  physically  important  noise  component  is  isotropic  noise  in  which  the  spatial 
density  of  the  noise  from  all  directions  is  the  same.  Since  there  is  no  preferred  direction, 
the  elements  of  the  isotropic  noise  cross -spectral  density  matrix  are  all  real.  In  a 
spherically  isotropic  field  the  cross -spectral  density  between  two  sensors  separated  by  a 
distance  d  is  sin(u)d/c)^(ud/c)  where  c  is  the  velocity  of  propagation.  The  symbol  Qq 
will  be  used  to  denote  this  isotropic  noise  cross -spectral  density  matrix.  Notice  that  Qq 
depends  only  on  the  geometry  of  the  array  and  frequency.  Another  isotropic  noise  field 
[14]  is  JQ(iod/c)  which  arises  in  two  dimensional  problems  when  the  field  is  cylindrically 
isotropic. 

Point  source  noise  comes  from  a  single  point  source  and  hence  is  signal-like  in  its 
spatial  characteristics. 

The  symbol 

Qd  =d(u)  d*(u)  (11) 

will  be  used  to  denote  this  type  of  noise  component.  Again  the  normalization  d*d  =  M  is 
assumed.  Isotropic  noise  is  simply  the  limit  of  summing  infinitesimal  independent  noise 
components  of  equal  intensity  from  all  directions. 

Henceforth,  in  order  to  simplify  the  notation,  the  dependence  of  various  quantities 
on  frequency  u>  will  usually  not  be  shown  explicitly.  * 

When  the  signal  matrix  has  the  simple  form  of(10)  the  expression  for  the  array  gain 
becomes 

~  k*Qk 


*  Many  of  the  results  expressed  in  general  vector-matrix  form  may  also  be  used  directly 
in  time-domain  formulation  of  the  problem  [7]  . 
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The  numerator  of  (12)  is  the  magnitude  squared  of  the  inner  product  of  the  steering  vector 
k  and  the  signal  direction  vector  m.  We  may  think  geometrically  in  terms  of  a  generalized 
angle  y  between  these  vectors  defined  as  follows: 


cos  y  = 


lr»r 


(k*k)  (m*m) 


Substituting  from  (13)  into  (8)  and  using  the  relationship  m*m  =  M  yields 

G  =  M  cos2  y  j  k*k  j(k*  Q  k) }  (14) 

The  array  gain  against  spatially  white  noise  is  therefore 

G  =  M  cos2y  (15) 

w 

which  is  also  the  normalized  signal  response  (numerator  of  (8)  ).  It  is  a  function  of  the 

misalignment  of  the  steering  vector  k  and  the  signal  direction  vector  m. 

In  a  conventional  beamformer  k  is  chosen  to  be  proportional  to  m.  This  choice  of 

"matching  to  the  signal  direction"  makes  y  equal  to  zero.  This  maximizes  the  gain  against 

* 

spatially  white  noise,  alias  normalized  signal  response.  The  maximum  value  of  G  in  (15) 

w 

is  M.  If  G^  is  less  than  unity  the  array  becomes  more  sensitive  than  a  single  sensor  to 
spatially  white  noise. 

The  factor  {  k*k  ^(k*  Q  k)  f  is  the  reciprocal  of  the  normalized  noise  response  (the 

denominator  of  (8)).  It  is  also  the  gain  ratio  G/G  which  we  shall  denote  by  p .  When  the 

w 

noise  is  spatially  white  this  factor  is  unity.  For  an  arbitrary  noise  matrix  this  gain  ratio 
p  is  bounded  as  follows: 

1  1 

-r -  <  -  -  —  (16) 

Xmax  k*  Q  k  Xmin 


where  X  .  and  X  are  respectively  the  smallest  and  largest  eigenvalues  of  Q.  Since 
min  max  — 

Q  is  normalized  to  have  its  trace  equal  to  M,  the  average  size  of  the  eigenvalues  of  Q  is 
unity. 

When  the  noise  is  isotropic  the  array  gain 

Gq  =Mcos2y  |  k.#k^(k  QQk)  (  (17) 

is  known  as  the  directivity  index.  The  gain  ratio 

G/G  =k*k/(k*Qk)  (18) 

is  known  as  the  Q-factor  T15,  24]  or  super-gain  ratio  f22] .  We  shall  call  p  the  general¬ 
ized  super  gain  ratio  since  it  involves  replacing  the  isotropic  noise  matrix  in  (18)  with  a 
general  noise  matrix.  The  isotropic  noise  covariance  matrix  depends  on  the  array  geome¬ 
try.  For  line  arrays  it  becomes  the  identity  matrix  when  elements  are  spaced  at  one-half 
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wavelength  intervals.  When  sensors  are  moved  closer  together  than  one-half  wavelength 
some  of  flie  eigenvalues  of  the  isotropic  matrix  become  extremely  small  and  p  may  be¬ 
come  large. 

"Optimum"  and  "adaptive"  beamformers  choose  k  to  be  something  other  than  pro¬ 
portional  to  m.  The  underlying  philosophy  is  to  accept  a  reduction  in  normalized  signal 
response  in  order  to  achieve  an  even  greater  reduction  in  normalized  noise  response, 
thereby  improving  the  gain.  Since  normalized  signal  response  and  gain  against  spatially 
white  noise  are  the  same,  these  processors  can  become  very  sensitive  to  white  noise  if 
too  great  a  reduction  in  normalized  signal  response  takes  place.  An  "optimum"  steering 
vector  k*  under  a  number  of  criteria  [7]  is 

k*  =  m*Q  _1  (19) 

Lewis  and  Schultheiss  [14]  present  an  interesting  discussion  of  the  gain  of  optimum  and 
conventional  processors  and  their  relationship  to  the  eigenvalues  of  the  noise  covariance 
matrix.  For  die  choice  of  k  given  by  ( 19)  the  gain,  white  noise  gain  and  generalized 


super -gain  ratio  become  respectively 
G  =  m*Q  *m 

(m*Q  -1  m)2 

°W  12*2  2  ID 


(20) 

(21) 


(22) 


The  ratio  (22)  of  course  still  satisfies  the  bound  given  by  ( 16).  McDonough  [  16] 
presents  a  complicated  expression  relating  (21)  to  the  ratio  of  the  largest  and  smallest 
eigenvalues  of  Q.  He  concludes  (21)  may  become  very  small  if  A.  /  X  .is  large. 

The  theory  of  super-directive  arrays  involves  achieving  a  high  directivity  index  by 
using  the  optimum  steering  vector  for  isotropic  noise  k*  =  m*^1.  The  literature  on 
super-directive  arrays  is  extensive  [6, 12, 22, 23, 24,25] .  It  is  well  known  that  high  direc¬ 
tivity  index  (nearly  equal  to  M^)  can  be  achieved  with  endfire  steering  of  arrays  of  closely 
spaced  elements.  In  order  to  achieve  these  high  gains  the  steering  vector  is  made  almost 
orthogonal  to  the  signal  direction  vector  and  the  gain  against  white  noise  becomes  ex¬ 
tremely  small.  For  example,  Gilbert  and  Morgan  [12]  report  a  maximum  directivity 
index  of  15.8  but  a  corresponding  white  noise  gain  of  only  1.5  x  10~4  with  a  four  element 
endfire  array  with  1/16  wavelength  inter-element  spacing.  They  also  present  the  result 
that  the  maximum  directivity  index  of  an  array  averaged  over  all  steering  directions  is 
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equal  to  the  number  of  elements  M.  To  see  this,  we  write  (20)  for  isotropic  noise  as 

G  =  trace  (Q"1  mm*)  (23) 

and  notice  that  the  average  of  mm*  over  all  directions  is  by  definition  the  isotropic  matrix 


Sensitivity  analyses  of  super -directive  arrays  have  found  the  super-gain  ratio  and 
the  reciprocal  of  the  white  noise  gain  to  be  useful  measures  of  sensitivity. 

In  this  paper  we  shall  find  that  1/G  and  p  play  key  roles  in  sensitivity  whenever 
the  perturbations  are  statistically  independent  from  element  to  element. 

The  causes  of  variations  of  P,  Q  and  k  from  their  nominal  values  can  be  many  and 
varied.  They  include:  violations  of  a  priori  assumptions  such  as  those  concerning  plane 
wave  or  point  source  nature  of  signal  or  interference,  or  isotropic  nature  of  the  noise 
field;  imperfect  measuring  devices  and  finite  measuring  intervals  used  to  estimate  Q;  non- 
stationarity  of  the  noise  field;  presence  of  signal  in  noise  estimates  in  passive  systems; 
position  errors  in  sensor  locations;  amplitude  and  phase  errors  in  analog  components;  and 
sampling  and  quantization  errors  in  digital  components.  In  adaptive  beamformers  the 
nature  of  these  variations  is  usually  related  to  the  inter -play  between  a  priori  assumptions 
and  on-line  measurements.  Sensitivity  to  variations  in  k  are  especially  important  since 
the  underlying  philosophy  of  adaptive  beamforming  is  continuing  adjustment  of  the  steering 
vector. 

m.  GENERAL  SENSITIVITY  RESULTS 
1.  Noise  Perturbations 

Suppose  that  the  noise  field  is  made  up  of  several  components  so  that  Q  maybe  ex¬ 
pressed  as  a  sum  of  Hermitian  matrices  of  the  form 

9  =  ai9i+<*222  + . +  an2n  (24) 

where 

n 

>  0,  i  =  1,  ••• ,  n  i  =  iai  =  1 

Each  Q.  in  (24)  is  non-negative  definite  and  has  its  trace  equal  to  M  so  that  it  is  a  legiti¬ 
mate  normalized  noise  cross-spectral  matrix.  The  parameter  a.  is  the  relative  strength 
of  the  i-th  component.  The  decomposition  of  Q  into  components  is  not  unique  but  is  some¬ 
times  useful  when  various  physical  sources  of  noise  exist.  When  Q  is  given  in  the  form  of 
(24)  the  array  gain  (7)  may  be  expressed  as  follows: 

G=(il1«i/Gi)"1  (25) 
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where  is  the  array  gain  against  the  i-th  noise  component;  that  is 
k*  Pk 


G.  = 


(26) 


1  k*  Qjk 

From  (25)  it  is  evident  that  a  large  value  of  a./G.  for  any  component  will  cause  the  gain  G 
to  be  small.  Thus  a  relatively  weak  noise  component  (small  &.)  can  limit  the  overall  array 
gain  if  the  gain  against  that  particular  component  is  very,  small.  Hence  choices  of  k  which 
give  low  gain  against  a  type  of  noise  \viiich  is  likely  to  be  present  result  in  systems  which 
are  sensitive  to  small  amounts  of  that  type  of  noise. 


Similarly  (6)  may  be  rewritten  as  follows: 

ct2 

Z  =  (J^  k*  Pk  [  1  H - j; 

s  — -  2. 

CT 

S 


.  £  a./G.  ] 

i  =  1  li 


(27) 


Quantitative  measures  of  sensitivity  to  changes  in  the  noise  matrix  can  be  obtained 
by  differentiating  (25)  and  (27).  Suppose  that 


Q  =  (1  -a)Q1+  <*Q2 


(28) 


The  fractional  sensitivity  of  G  and  z  to  the  substitution  of  a  little  Q2 -type  noise  for  an  equal 


amount  of  Q^-type  noise  are 


=1-«VG2> 


and 


or 


dz/da  _  _ 

'  z  a  - o  1  +  G 


(Gj/G^  -  1 


(29) 


(30) 


.dG/da. 

dz/da  _  G  oc=  o 
(  z  }a  =  o~  1  +  (S/N) 


(31) 


2  2 

where  (S/N)  =(0/0  )G,  is  the  output  signal  to  noise  ratio.  WhenG./G_  >  1,  the  sub - 
o  s  n  1  1  * 

stitution  of  Q2 -type  noise  for  Q^type  noise  causes  a  decrease  in  gain  and  an  increase  in 
output  power.  The  gain  ratio  in  (29)  and  (30)  may  be  expressed  as 


Gi  **92!l 
°2 =  srSiiL 


(32) 


In  the  special  case  of  spatially  white  noise  (32)  reduces  to  the  generalized  super-gain  ratio 
G, 


k*k 

~  — -  =  P 

k*  Q,  k 


(33) 


w  *  yji 


8 


MPL-U  -57/72 


discussed  earlier.  If,  in  addition  Q^-type  noise  is  isotropic,  (33)  becomes  the  super -gain 
ratio  or  Q- factor. 


2.  Signal  Perturbations 

Similarly  we  may  conceive  of  the  signal  field  being  made  up  of  several  components 


so  that  P  may  be  expressed  as 
P  = 


n 

V 


^  p.  P.  ;  p.  >  0  i  =  1,  *•  •  ,  n  ; 
i  =  l  l—i  l 


i=  i  h  =  1 


(34) 


Each  Pj  is  Hermitian  and  non-negative  definite  with  trace  equal  to  M. 


Then 


n 


and 


where 


G  =  (3.  G. 

x  =  l  l  i 


z  =  a  k*Qk  [l  +  — r- 

n  — -  ^ 

a 

n 


k*  P.  k 

G.  —  ■■ . ■■■ . 

1  k*  Q  k 


•  S  ,  P.G.] 

1  =  1  ri  l 


(35) 

(36) 

(37) 


is  the  gain  for  the  i-th  component  of  the  signal  field.  The  overall  gain  is  simply  the  sum 
of  the  gain  for  each  component  weighted  by  its  relative  strength. 

Suppose  that 

P  =(1  -  p  )P1+  p  P2  (38) 

The  fractional  sensitivities  of  G  and  z  to  the  substitution  of  a  little  ?2  -type  signal  for  an 
equal  amount  of  P^ -type  signal  are 


,  dG/d  p . 

K  G  p  =o 


-  1 


and 


|3-0 


<°2^1>  -  1 
i+  °^/<Gl  %> 


or 


,dz/dB  (S/N>o 

'  z  ;p=o  1  +  (S/N) 


( dG/d  p . 

{  G  '  p  =o 


(39) 


(40) 


(41) 


The  gain  ratio  G2/Gj  appearing  in  (39)  and  (40)  maybe  expressed  as 


°2  r  f2  - 

°1  k*  ?i  k 


(42) 
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When  is  the  result  of  a  point  signal  source  with  direction  vector  m  and^  *s  the  identity 
matrix  representing  independent  signal  perturbations  from  sensor  to  sensor  (42)  becomes 


VG1  =  1/GW 

(43) 

G2  =  p 

(44) 

which  are  by  now  familiar  expressions.  Low  white  noise  gain  in  (43)  means  large  response 
to  a  spatially  incoherent  signal  component. 

Before  proceeding  it  is  useful  to  consider  how  the  situation  represented  by  (38) 
might  arise.  Suppose  that  the  signal  vector  suffered  random  perturbations  so  that  it  could 
be  written  as 

m  =  \!  1  -  P  m^  +  s/p  6  (45) 

when  6  is  a  vector  of  random  perturbations  with  zero  mean  and  normalized  covariance 
matrix  E[6o*]  =  P^.  The  factor  \l  -(3  in  (45)  provides  the  normalization  so  that 

E  [m*m]  =  trace  [  E(mm*)]  =  M  (46) 


Then 

E[|k*m|2]  =  (1  -  P )  |  £*2^  | 2  +  p  k*P2k  (47) 

In  this  situation,  the  sensitivity  results  (39)  and  (40)  should  be  interpreted  on  an  ensemble 
average  basis. 

From  (39)  it  is  evident  that  small  perturbations  on  the  average  will  not  cause 
severe  degradation  in  G  or  z  even  if  «  G^.  However,  when  G^  »  Gj  rapid  increases 
in  G  and  z  are  to  be  expected.  Of  course,  any  individual  perturbation  of  the  type  (45)  may 
deviate  significantly  from  this  average  behavior.  From  (43)  we  see  that  when  k  is  nearly 
orthogonal  to  m  so  that  1/GW  is  large,  random  signal  perturbations  will  on  the  average 
cause  an  increase  in  G  and  z  by  decreasing  y .  Of  course,  unpredictable  increases  in  G 
and  z  due  to  slight  perturbations  in  m  with  no  change  in  signal  power  can  be  a  source  of 
confusion  and  may  be  undesirable.  This  is  especially  true  in  processors  which  are  sup¬ 
posedly  constrained. 

Another  approach  to  the  formulation  of  the  signal  sensitivity  problem  is  to  assume 
P  is  of  the  dyadic  form  of  (10)  and  to  allow  amplitude  and  phase  perturbations  to  occur  to 
the  individual  components  of  m,  i.  e. , 


m.  =  m?(l  +  a.)  exp  (\M  £.)  (48) 

)  J  ]  J 


turbations  have  zero  mean  and  are  assumed  to  be  independent  of  each  other.  This  approach 
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has  been  taken  in  a  recent  paper  by  McDonough  [16]  following  the  lead  of  Gilbert  and  Morgan 
[12] .  When  m.  is  given  by  (48),  the  quantity  E(  |  k*m  |  becomes 

3  .  10  «  rv 


E(  I  k*m  I  2)  =  k?k.m°m°*  C.. 

I - '  i,j  1  j  1  j  i] 


where 


'  1+E  (ah. 


C..  = 
i] 


for  i  =  j 


(49) 


(50) 


|  1  +  E(a.a^)  |  E[exp(  \ -1  (£.  -  £  ))],  for  i  /  j 

Notice  that  the  quantity  m  in  (48)  is  not  normalized  as  m  was  in  (46)  since  for  m  defined 
in  (48) 

E  [m*m]  =M+  Z  |m.|^  E(a^)  >  M  (51) 

-  j  J  J 

As  shown  by  McDonough,  (49)  may  be  simplified  under  the  following  additional 
assumptions  which  were  made  by  Gilbert  and  Morgan: 

1.  Amplitude  perturbations  are  independent  from  sensor  to  sensor  and  of  equal 

variance,  i.e. ,  E[a.  a.]  =  6..  . 

1  1  a  ij 

2.  Phase  perturbations  are  Gaussian,  small,  independent  from  sensor  to  sensor 

2  2 

and  of  equal  variance,  i.e.,  E[£.|.]  =  ° {.  <5..,  ° ^  «  1. 

ty  3  x3 

3.  m.  is  a  pure  phasing,  i.e. ,  |  m.  |  -  1 ,  j  =  1 ,  •  •  •  »  n. 

Under  the  above  assumptions 


C.. 

13 


=  1+  <T 


=  1-0, 


for  i  =  j 
for  i  /  j 


and  (49 )  becomes 


E[  |  k*m  |  ^  ]  =  (1  -  o^  )  j  k*m°  j  ^  +  (o^  +  0^  )  k*k 


(52) 


(53) 


which  is  almost  identical  to  (47)  for  the  corresponding  case  £2  =  I_.  The  difference  can  be 
attributed  to  the  normalization  of  (45). 

In  general  (49)  can  be  made  to  look  like  (47)  by  defining  a  matrix  F  with  elements 


F..  =  m.  mT  (C..  -  1) 
13  1  3  13 


Then  (49)  may  be  written 


E  [  j  k*m  | 1  ]  =  |  k*mj  2  +  k*Fk 


(52) 


(53) 


Much  of  the  sensitivity  work  [12, 16, 24]  in  the  field  of  antenna  arrays  has  been  con¬ 
cerned  with  the  quantity  E T  |  k"m  |  ^  ]/|  k^ mQ  |  ^  when  the  components  of  k  and/or  m  have 
been  perturbed  in  amplitude  and/or  phase  from  their  nominal  values  kQ  and  mQ.  From 
(49)  it  is  evident  that  the  type  of  result  will  be  the  same  whether  C„  arises  from 
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perturbations  of  m,  k,  or  some  combination  of  both. 


3.  Steering  Perturbations 

As  suggested  above  there  are  similarities  between  signal  perturbations  and  per¬ 
turbations  in  the  steering  vector  k.  Perturbations  in  k  may  also  be  approached  in  two  ways, 
analogous  to  (45)  and  (48)  respectively.  However,  there  is  an  important  difference  in  that 
perturbations  in  k  affect  both  signal  and  noise  terms. 

Following  the  approach  of  (45),  suppose  that 

k  =  \'l  -  (  k  1  +  \'T  77  (54) 

where  t)  is  a  vector  of  random  perturbations  with  zero  mean  and  covariance  matrix 
E[t}77*]  =  W  where  W  is  normalized  to  have  its  trace  equal  to  k*kj.  Then  expected  out¬ 
put  power  is 

z  =  E[k*  Rk ]  =  (1  -  e  )k*Rk  +  f  trace (RW)  (55) 

and  the  sensitivity  of  the  expected  output  power  is 
Av./Af  trace  (RW) 


.dz/de 

{  z  6  =  0 


k  jR  k  L 


Defining  a  gain  G  as  follows: 

G  =  E(k*Pk)/E(k*  Q  k) 


we  obtain 


(1  -  e)k*  Pk1  +  e  trace  (PW) 


(1  -  e)  k*  Q  k  1  +  e  trace  (Q  W) 

The  sensitivity  of  this  gain  is 

»  trace  (PW/Gj)  -  trace  (QW) 


,  dG/d  e 

1  5  €=0 


!i*l2  111 


Simplifications  again  occur  when  the  perturbations  are  independent  from  sensor  to  sensor 
and  of  equal  variance  so  that  W  =  (k*  k  /M  )^  •  Then  (55)  becomes 

_  _  /I  \  1  i  n  1  ■  _  i _ 2  _ 2  .  ,  *  , 


z  =  (1  -  e )  k*  R  k  ,  +  e  (ct  +  a  )  k*  k  . 

-1 - 1  s  n  -1-1 


and  (56)  becomes 


. dz/df 
'  z  'e  =0 


k*  k 
5rl  Ll 

k*l  Q  k  1 


1  +  <S/N\n 
1  +  (S/N). 


When  (S/N).^  and  (S/N)q  are  both  much  less  than  unity  the  output  power  sensitivity  of  (61) 
is  approximately  equal  to  p  -  1.  When  (S/N).n  and  (S/N)q  are  both  much  larger  than  unity,  it 
is  approximately  equal  to  |  (1/Gw)  -  1  |  . 

A  similar  simplification  occurs  in  (58)  which  becomes 
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G  = 


(1  -  c)k*  Pkj  +  f  k*  kl 
<1-f>k*1Qki+ 


and  (59)  which  becomes 

=3 


__  Jr*  Jr 

dG/de  _  -1  -1 
(  G  >€=0  k*Qk 


1 

G. 


-  1 


(62) 


(63) 


—  i  —  — 1  <  1 

For  the  usual  case  of  G^»  1  the  gain  sensitivity  of  (63)  is  approximately  equal  to  -  P  . 

Following  the  alternate  approach  of  (48)  and  considering  amplitude  and  phase  per¬ 
turbations,  we  may  describe  k  as  follows: 

kj  =  k°  (1  +  bj)  exp  (\^1  0^)  (64) 

where  k°  is  the  nominal  value  of  the  j  -th  component  of  k.  The  amplitude  perturbation  b^ 

and  the  phase  perturbation  <p.  have  zero  mean  and  are  assumed  to  be  independent  of  each 

other.  When  k  is  described  by  (62)  the  expected  output  power  is 

z  =  E(k*Rk)  =  I  kf*  k?  D..  R.. 

-  l,j  1  )  »  'J 

where  is  the  i,  j-th  component  of  R  and  in  analogy  to  (50) 


(65) 


D«  = 


1  +  E(b7)  , 


for  i  =  j 


1 1  +  E(b.b.)  |  E[exp  (  'sTT  (0.  -  0. )  )  ]  ,  for  i  /  j 


(66) 


Similarly  G  defined  in  (57)  becomes 


G  =  ill 


Z  k?*  k?  D..  P„ 


j  -ij  ‘ij 


Z  k°*k°D..Q.. 
i  J  i]  ij 


(67) 


M 


where  P  and  Q_  are  i,  j-th  components  of  P  and  Q  respectively.  Under  the  following 
assumptions  (63)  and  (65)  maybe  simplified  considerably: 

1.  Amplitude  perturbations  are  independent  from  sensor  to  sensor  and  of  equal 

2 

variance,  i.e.,  E[b,b.]  =  ct,  6..  . 

i  J  b  ij 

2.  Phase  perturbations  are  Gaussian,  small,  independent  from  sensor  to  sensor 

2  2 

and  of  equal  variance,  i.e.,  E[0  0  J  =  Oa  6..  ,  Ga  «  1  . 

12  I  1 2  1  J  9  l)  9 

=  |  k.  j  for  i,  j  =  1,  •  •  •  M,  or  =  P^  =  1  for  i  =  1,  •  •  • ,  M. 


Then  (65)  becomes 


z  =  i°*M°  +  <CT0  +  C7bHCTs  +  CTn)!iO*KO 


(69) 


which  closely  resembles  (60).  Under  the  same  assumptions  (67)  becomes 
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G  = 


(1  -  o2^)  k°*Pk0  +  (a^  +  <J2)  k°+  k° 
(1  -  k°*  Qk°  +  (ct^  +  02)  k°*  k° 


(70) 


which  closely  resembles  (62). 

Finally  let  us  note  that  when  m  is  perturbed  as  in  (48)  and  k  is  independently  per¬ 
turbed  as  in  (64), 


El 

k*  m  2] 

E 

Ik*  Q  k] 

£  k°*  k°  m°  m°*  D..  C.. 

i  ]  i  J  U  i] 

£  k°*  k°  D..  Q.. 

i  J  iJ  ij 


(71) 


4.  Summary 

Table  I  presents  a  summary  of  some  of  the  more  important  sensitivity  results 
derived  in  this  Section. 


IV.  INTERFERENCE  REJECTION  AND  SIGNAL  SUPPRESSION 

The  problem  of  designing  an  array  to  be  insensitive  to  a  point  source  of  noise  has 
received  considerable  attention  in  the  literature  [l,  2,  8, 20] .  In  the  notation  of  this  paper, 
the  goal  is  to  reject  a  component  of  the  noise  field  of  the  form  of  (11). 

The  gain  of  any  beamformer  against  such  an  interference  is 

=  |  k*m|  |  k*l| 2  =  cos2y  jcos2Q  (72) 

where  0  is  the  generalized  angle  between  the  steering  vector  k  and  the  interference  direc- 

2 

tion  vector  d.  The  quantity  cos  0  is  obtained  by  replacing  m  in  (13)  by  d. 

Perhaps,  the  conceptually  simplest  approach  to  interference  rejection  is  to  use  a 

2 

steering  vector  which  is  orthogonal  to  d  so  that  cos  0  =0.  Such  a  beamformer  com¬ 
pletely  nulls  out  the  unwanted  interference.  The  quality  of  the  null  will  be  degraded  if 
either  k  or  d  suffers  random  perturbations.  Because  of  the  symmetry  of  the  expression 

cos20  =  |k*dj2^  (k*k)  (d*d)  (73) 

perturbations  in  k  or  d  have  similar  effects  on  the  null.  Perturbations  in  k  also  affect  the 
signal  response  while  perturbations  in  d  do  not.  The  quantity  E(  |  k*d  j 2),  when  d  suffers 
random  perturbations,  maybe  obtained  directly  from  the  results  (47)  or  (49)  presented  for 
the  E(  J  k*m  j2),  when  m  underwent  random  perturbations.  In  the  simplest  case  of  inde¬ 
pendent  perturbations  from  sensor  to  sensor 

E[cos20  ]  =  a2/ M  (74) 
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where  cr  is  equal  to  |5  in  (47)  or  the  combined  variance  of  the  amplitude  and  phase  per¬ 
turbations  in  (53).  In  this  situation  the  expected  gain  is 

15(G  )  =  G  /a2  (75) 

d  w 

Thus  ifG^  is  small,  nulling  becomes  impractical. 

'Hie  null-steering  beamformer  of  Anderson  [l,2]  uses  the  following  steering  vector  " 
k*  =  m*  [i  -dd'/M]  (76) 

For  this  steering  vector  the  white  noise  gain  G  is  given  by  the  following  equation: 

G  =  M[l  -  jm  dJ^/M^]  =  M[l  -  cos^g  ]  =  M  sin^g  (77) 

where  g  is  the  generalized  angle  between  the  signal  direction  vector  and  the  interference 
direction  vector.  Thus,  the  expected  gain  (75)  against  the  interference  will  be  low  when  the 
interference  and  signal  are  closely  spaced  such  as  when  the  interference  is  within  the  main 
lobe  of  a  beam  steered  conventionally  in  the  signal  direction. 

When  an  interference  is  added  to  an  existing  noise  field  the  noise  matrix  becomes 


Q  =  CTi  Qj  +  11*  =  <CT1  + 


The  optimum  steering  vector  is 


Then 


and 


k*  =  m* Q1  =  [1  +  <ct^/ct^)]  m* 


Q.i  siVs;1 


k*d  =  [1  +  (ct^/o^)]  m*  Q^d 


k*m  =  [l  +  fcrj/o2)]  m* Q^1  m 


1  - 


-1  2  2 
l*Qj  I+Oj/cTj 


i  +  [d*Q11  d  a^/(7j]  1 


1  - 


Ibtsi1!!2  1/ _ 1 _ \) 

.  \ 1  +  1  ad/cri 1  1  /  i 


(78) 


(79) 

(80) 


(81) 


The  gain  of  the  "optimum  beamformer  is  also  given  by  (81). 

In  order  to  interpret  (80)  and  (81)  it  is  useful  to  consider  the  meaning  of  various 
quantities  which  appear  in  these  equations: 

2  2 

a  ,/a .  --  interference  to  noise  ratio  at  the  input  (I/N).  . 
d  1  in 

*  The  matrix  [j.  -  dd*/M]  is  a  projection  operator  which  passes  only  the  component  of  m 
which  is  orthogonal  to  d. 
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m  Qj  =  ^  --  optimum  steering  vector  in  the  absence  of  interference. 

-1  2,2 

d  Q,  d  a~/a,  =  (1/N)  --  maximum  possible  output  interfererce-to -noise 

—  —  l  —  d  1  max 

ratio,  that  is,  the  output  interference-to -noise 
ratio  of  an  optimum  beamformer,  for  -type 
noise,  steered  in  die  direction  of  the  inter- 

!!  ference. 

1 2 

k^d  —  response  in  die  interference  direction  of  an  opti¬ 
mum  beamformer,  for  -type  noise,  steered 
in  the  signal  direction. 

m  ■"  m  --  optimum  gain  in  the  absence  of  interference. 

The  ''optimum”  beamformer  does  a  tradeoff  between  nulling  the  interference  and 

preserving  gain  against  -type  noise.  This  tradeoff  involves  the  interference-to-noise 

ratio  and  the  sidelobe  level  in  the  interference  direction  of  the  optimum  beamformer,  for 

Q^-type  noise,  steered  in  the  signal  direction.  When  Q^-type  noise  is  spatially  white  so 

that  Qj  M,  the  limit  of  the  steering  vector  (79)  for  large  input  interference-to-noise  ratios 

is  die  Anderson  null-steering  processor  (76)  to  within  a  constant  of  proportionality. 

When  (1/N)  is  small  the  interference  has  little  effect  on  the  optimum  steering 

vector  which  remains  nearly  the  same  as  it  would  be  in  the  absence  of  interference.  When 

(I/N)  becomes  large,  k*d  given  by  (80)  becomes  small  and  a  null  develops  in  the  direc- 
max  ■ 

tion  c  f  the  interference.  The  factor 


appearing  in  (81)  must  be  less  than  or  equal  to  unity  by  the  Schwarz  inequality.  It  may  be 
interpreted  as  the  cosine  squared  of  the  generalized  angle  between  m  and  d  in  the  linear 
vector  space  in  which  length  is  defined  relative  to  the  metric  so  that  is  the 

length  of  m. 

The  question  of  signal  suppression  arises  when  the  steering  vector  k*  is  made  pro¬ 
portional  to  m*  R  *  instead  of  m*Q  1  and  R  1  contains  the  signal  vector.  It  is  of  particu¬ 
lar  concern  since  the  nulling  embodied  in  the  R  *  operation  usually  involves  measured 
quantities  and  die  matching  operation  involved  in  m*  usually  is  based  on  a  priori  assump¬ 
tions  about  the  signal  direction  vector. 

Mathematically,  signal  suppression  may  be  treated  as  a  special  case  of  interference 
rejection  so  that  the  results  developed  above  may  be  applied  directly.  In  particular,  we 
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may  look  at  the  effect  of  the  inclusion  of  signal  in  the  estimator  [4] 

z  =  ("fR-1™)’1  (83) 

which  is  the  power  output  of  a  beamformer  with  k*  =  m*  R  ^(m*R  *m).  This  k  provides 

a  minimum  variance  unbiased  linear  estimate  of  the  signal  [4,7].  The  quantity  z  ^  = 

-1  2  2  2 
m"  R  m  may  be  obtained  by  replacing  the  initial  factor  (1  +  cr^/a  ^ )  in  (81)  with  (l/o^) 

and  redefining  the  other  terms  in  (81)  as  follows: 

d:  actual  (measured)  signal  direction  vector 

m:  assumed  (a  priori)  signal  direction  vector 

Q^:  noise  only  matrix 

Oy  input  noise  level 

o^:  input  signal  level 

Mismatch  occurs  when  m  /  d.  The  effect  of  mismatch  can  be  seen  by  examining  the  ratio 


(2Wd 

rQl1! 

<Z)m  =  d 

m^Qj1 

-1  2  2-1 
l  +  td’Qj  d  cd/o\) 


‘Q;1>  (l+rs/d  tyy 


There  are  two  distinct  effects  of  mismatch:  First  is  the  effect  of  the  factor 

This  effect  is  the  usual  effect  of  mismatch  discussed  in  Section  III 

and  has  nothing  to  do  with  the  inclusion  of  the  signal  in  the  matrix  inversion.  Second  is 

2 

what  we  shall  call  the  anomalous  signal  suppression  caused  by  die  presence  of  cos  (Mn_i) 

^1 

in  (84)  which  is  a  direct  result  of  including  the  signal  in  the  matrix  inversion  process. 

The  anomalous  signal  suppression  will  be  insignificant  as  long  as  the  quantity 
-1  2  2 

(cP'Qj  f?  CTd^CTl  ^  *S  sma^>  as  l°ng  as  die  signal -to -noise  ratio  at  the  output  of  a 

perfectly  matched  optimum  beamformer  with  k*  =  d *  Q ^  would  be  small.  If  this  weak  sig¬ 
nal  criterion  is  not  satisfied,  the  ratio  given  in  (82)  will  be  reduced  through  the  inclusion  of 
the  signal  in  the  matrix  inversion  and  the  mismatch  between  the  actual  and  assumed  signal 
direction  vectors. 


-1  2  2 

If  we  define  (S/N)max  as  (d*Q^  d  a  ^/a  ^  )  then  (84)  maybe  simplified  to  the 


following: 
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(7Wd  dQ^d 

(z>m=d 


1 


1  +  (S/N)max  Sin  VlJ 
1 


(«5) 


2  2  2  -1 
where  sin  (• )  =  1  -  cos  (• ).  The  factor  [l  +  (S/N)  sin  (p_  i  )]  is  a  direct  measure 

max  1 

of  the  anomalous  signal  suppression. 


V.  OPTIMIZATION 

There  are  a  number  of  closely  related  optimization  problems  which  have  been  or 
can  be  formulated.  In  this  section  we  shall  briefly  sketch  some  of  these  problems  with 
emphasis  on  formulations  which  somehow  take  sensitivity'  into  consideration. 

A.  Beam  formers 

1.  Maximum  Gain  and  Minimum  Variance 

The  problem  of  unconstrained  array  gain  maximization  is  to  choose  k  such  that 
G=k*Pk/k*Qk  (7) 

is  maximized.  The  solution  of  this  problem  is  known  [6]  to  be  choosing  k  proportional  to 
the  eigenvector  corresponding  to  the  largest  eigenvalue  of  (Q  1 P).  When  r  =  mm*,  the 
optimum  k  is  proportional  to  Qm  which  is  well  known  [7,  10,  17]  and  easily  shown  by  direct 
application  of  the  Schwarz  inequality  [9] . 

The  relationship  between  maximizing  the  gain  and  minimizing  variance  under  a  con¬ 
straint  on  signal  response  lies  in  that  maximizing  G  is  equivalent  to  minimizing  the  de¬ 
nominator  of  (7)  subject  of  a  constraint  on  the  numerator.  Since  the  gain  in  (7)  is  invariant 
to  a  scaling  of  k,  the  imposition  of  a  constraint  of  the  form  k*m  =  1  has  no  effect  on  G  and 
simply  determines  the  constant  of  proportionality.  Thus  k  =  Q  2}/(nj*Q  Jm)  both  maxi¬ 
mizes  G  and  satisfies  the  constraint  k*m  =  1.  The  related  problem  of  minimizing  k*Rk 
subject  to  the  same  constraint  leads  to  k  =  R  Srn  / (m*  R  *m). 

2.  Maximum  Gain  Subject  to  Sensitivity  Constraint 

As  shown  in  Section  III  the  generalized  super-gain  ratio  p  and  the  gain  against 
white  noise  G^  play  key  roles  in  the  sensitivity  of  beamformers  when  the  perturbations  are 
independent  from  sensor  to  sensor.  The  problem  of  maximizing  the  gain  (7)  subject  to  a 
constraint  on 

Gw  =  (86) 
is  most  easily  formulated  as  that  of  finding  the  k  which  minimizes  1/G  with  a  constraint  on 
1/Gw>  That  is  minimizing 
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1/G+  X/G  =  k  (Q  UI)k  /k^Pk  (87) 

where  \  is  a  Lagrange  multiplier.  Finding  the  k  which  minimizes  (87)  is  equivalent  to 
finding  the  k  which  maximizes  the  reciprocal  of  (87),  which  in  turn  is  an  eigenvalue  prob¬ 
lem  of  tile  type  discussed  above.  Hie  optimum  k  is  chosen  to  be  proportional  to  the  eigen¬ 
vector  corresponding  to  the  largest  eigenvalue  of  [Q  +  \  l]  .  This  problem  was  first 
addressed  by  Gilbert  and  Morgan  [12]  with  P  =  m  ni*  so  that  the  optimum  k  is  equal  to 
[Q  +  A  j]  . 

A  related  problem  is  to  maximize  G  subject  to  a  constraint  on  p  .  In  this  case  v/e 
maximize 

G  +  Ap  =  k*[P  +  AJ]  k/k  Q  k  (88) 

Again  we  have  a  similar  eigenvalue  problem.  Uzsoky  and  Solymar  [24]  examined  this 
problem  for  the  case  of  isotropic  noise  Q  =Qq.  For  the  formulation  of  (88)  to  be  consistent 
p  must  be  specified  within  its  allowable  range  given  by  (16).  From  (86)  we  see  that  the 
effect  of  adding  the  constraint  on  p  is  equivalent  to  that  of  adding  a  spatially  white  (spa¬ 
tially  incoherent)  component  to  the  signal  field. 

Lo,  Lee  and  Lee  1 15]  review  a  number  of  optimization  problems  and  present  a 
numerical  approach  to  the  solution  of  the  more  difficult  problem  of  maximizing  G  subject 
to  a  constraint  on  the  Q -factor.  They  maximize 

<  | | 2/  k*Qk)  +  <Ak’k/k*Qok) 

All  of  these  optimization  problems  are  insensitive  to  a  scaling  of  k.  Hence  the 
addition  of  a  linear  constraint  such  as  k^m  =  a  /  o  is  handled  by  simply  scaling  k. 

3.  Maximization  of  Expected  Quantities 

A  somewhat  different  approach  to  the  problem  of  sensitivity  is  to  take  the  type  and 
anticipated  size  of  perturbations  into  account  before  maximizing  G.  Thus  perturbation 
terms  are  included  in  the  ratio  to  be  maximized.  A  special  case  of  this  approach  has  been 
taken  by  Cheng  and  Tseng  [6]  who  maximize  (71)  and  present  numerical  results  for  a  linear 
endfire  array  of  eight  dipoles.  Again  the  magnitude  of  k  remains  free  so  that  an  additional 
constraint  of  the  form  k'm  =  a  ^  o  can  bn  handled  by  scaling  k. 

4.  Multiple  Linear  Constraints 

In  the  preceding  discussion  a  single  linear  constraint  could  usually  be  handled  by 
using  the  degree  of  freedom  of  the  magnitude  of  k  which  was  left  undetermined  in  the  pro¬ 
cess  of  gain  maximization.  An  exception  arises  when  the  constraint  is  a  null  of  the  form 
k*d  =  o.  This  situation  may  be  treated  as  a  special  case  of  the  more  general  problem  of 
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multiple  linear  constraints.  Multiple  constraints  may  be  used  to  reduce  sensitivity  to  sig¬ 
nal  perturbations  by  keeping  the  gain  relatively  constant  over  a  range  of  signal  perturba¬ 
tions.  A  similar  technique  may  be  used  to  control  sidelobes  in  a  specific  neighborhood  of 
directions.  The  problem  can  be  formulated  as  that  of  minimizing  the  output  power  z  subject 
to  the  constraint  H*k  =  g.  Each  row  vector  of  the  constraint  matrix  H*  imposes  a  con¬ 
straint  of  die  form  h*k  =  g..  Thus,  the  constraint  matrix  H*  has  a  row  for  each  constraint. 
— x  —  i  - 

The  total  number  of  rows  must  be  less  than  the  number  of  sensors  or  die  problem  will  be 
overspecified.  Using  a  Lagrange  multiplier  vector  A*  we  may  minimize 

z  =  k*Rk  +  X*[H*k  -  g]  +  [k*H-g*U  (89) 

Completing  the  square  [7]  yields 

z  =  [k*  +  A*  H*  R_1]  R[R-1H^+k]  -  X*H*R-1H\  -\*g-g*A  (90) 

Since  k  appears  only  in  the  initial  quadratic  term  of  (90),  die  solution  is  obviously  to  make 
that  term  equal  to  zero  by  choosing 

k  =  -R_1HA.  (91) 

Using  the  constraint  H*k  =  g  to  eliminate  A  ,  finally  yields 

k  =  R-1H  [H*  R'1  H]"1  g  (92) 

For  this  choice  of  k  the  output  power  z  becomes 

z  =  g*  [H*  R-1h]  ”^g  (93) 

The  use  of  an  estimate  for  R"1  in  (91)  is  a  generalization  of  the  estimator  (83).  One  ap¬ 
proach  to  maintaining  signal  response  is  to  have  each  row  vector  of  H*  be  a  steering  vector 
in  the  neighborhood  of  the  steering  direction  m  and  to  let  g  be  a  vector  with  the  number  one 
as  each  component. 

B.  General  Array  Processor 

So  far  our  discussion  has  centered  on  beamformers  which  have  the  structure  of 
Figure  1.  A  more  general  array  processor  structure  is  illustrated  in  Figure  2.  In  the 
more  general  processor  there  are  multiple  outputs  obtained  by  a  matrix  filtering  of  the 
input  signals.  In  the  processor  of  Figure  2,  K*  is  a  matrix.  While  the  processor  of  Fig¬ 
ure  1  could  be  followed  by  a  simple  square-law  detector  and  averager,  the  matrix  filter  K* 
may  be  followed  by  a  more  general  quadratic  processor.  It  is  known  [7,  18]  that  the 
structure  of  Figure  1  is  only  optimum  when  P  is  a  simple  dyadic.  Since  the  effects  of  per¬ 
turbation  is  to  destroy  the  dyadic  nature  of  the  expectation  of  P  we  are  naturally  led  to  the 
consideration  of  more  the  general  processor. 
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For  example,  suppose  that  the  signal  and  noise  were  both  completely  incoherent 
from  sensor  to  sensor,  i.e. ,  P  =  Q  = J.  Then  all  beamformers  of  the  type  of  Figure  1  will 
provide  no  gain  (G  =  1)  while  an  incoherent  combination  of  the  sensor  outputs  can  provide  a 
gain  of  \rM. 

A  processor  which  forms  multiple  closely  spaced  beams  and  averages  output  power 
across  these  beams  is  a  special  case  of  the  general  array  processor. 

In  the  general  array  processor  the  output  power  is  defined  as 

z  =  trace  (K*  R  K)  =  E(y)  (94) 

In  order  to  be  able  to  properly  account  for  the  potential  of  incoherent  gain,  we  define  a  gain 
in  terms  of  the  detection  index  at  the  output  of  a  general  quadratic  processor  [7]. 


G  = 


E  (y  i  S  +  N )  -  E(y|  N)  , 

t<\ 

|E(y2  |  N)  -  E2(y  1 N)  j 

|l/2  \ 

w2s/ 

(95) 


When  the  noise  is  Gaussian  [?],  (95)  may  be  written  as 

G  =  trace  [K*  PK  ]  /{ trace  [(K*  Q  K)2  ]  {  1/2  (96) 

Notice  that  (96)  reduces  to  (7)  when  K  is  a  column  vector. 

The  problem  of  maximizing  (96)  has  been  solved  [7,  9]  by  direct  application  of  the 
following  Schwarz  inequality: 

1 2 


j  trace  (A*  B)  J  <  trace  (A*  A)  trace  (B*  B) 
The  optimum  choice  of  K  is 


K  =  cQ'1  A 


(97) 


(98) 


where  c  is  an  arbitrary  scalar  constant  of  proportionality  and  P  =  AA*.  The  matrix  A, 
assumed  for  convenience  to  be  of  hill  rank,  will  have  M  rows  and  r  columns  where  r  is  the 
rank  of  the  matrix  P.  The  maximum  value  of  gain  (96)  is 

G  =  {trace  [(PQ"1)2]  }  1/2  (99) 

This  maximum  gain  cannot  be  achieved  by  any  K  with  less  than  r  columns. 

The  problem  of  maximizing  gain  (96)  determines  K  to  within  a  scalar  multiple. 

Thus,  constraints  of  the  form  trace  (H*  K)  =  a  ^  o  may  be  handled  by  a  simple  scaling  as 
before. 

The  problem  of  minimizing  trace  (K*  R  K)  subject  to  multiple  linear  constraints  of 
the  form  H*  K  =  L  is  readily  handled  by  completing  the  square  similar  to  the  way  it  was 
done  in  (89).  For  details,  see  the  development  in  equations  (28)  to  (31)  of  [7].  The  solution 
is 

K  =  R_1H  fH*  R"1  H  j  _1  L  (100) 
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for  this  value  of  K, 

K  *  R  K  =  L*  [if  R_1  H  ]  _1  L  (101) 

Equations  (100)  and  (101)  are  generalizations  of  (92)  and  (93).  A  further  generalization  of 
the  estimator  (S3)  is  therefore 

z  =  trace  \  L*  [H*  R-1  H]  _1  L  \  (102) 

C.  Implementation 

Most  of  the  work  [5, 13, 21, 26]  in  adaptive  beamforming  has  dealt  with  unconstrained 
optimization  problems  so  that  there  remains  much  unbroken  ground  in  the  field  of  con¬ 
strained  optimization. 

Some  exceptions  do  exist.  The  work  in  antenna  array  optimization  [6,  12,22,24]  has 
been  concerned  with  sensitivities  but  not  concerned  with  appropriate  algorithms  for  on-line 
adaptation.  The  problem  of  multiple  linear  constraints  leads  naturally  to  stochastic  ver¬ 
sions  of  gradient -projection  type  algorithms  [l  1 , 19] .  This  approach  has  also  been  sug¬ 
gested  [27]  for  non-linear  constraints.  However, slow  convergence  maybe  anticipated. 

VI.  CONCLUSION 

General  sensitivity  measures  have  been  developed  for  the  cases  of  perturbations  to 
the  signal  field,  the  noise  field  and  the  steering  vector.  These  sensitivity  measures  may 
be  used  to  test  the  practicality  of  particular  processors.  It  was  found  that  the  white  noise 
gain  and  the  generalized  super -gain  ratio  are  key  parameters  in  determining  the  various 
sensitivities  when  the  perturbations  are  independent  from  sensor  to  sensor. 

The  problem  of  anomalous  signal  suppression  through  the  inclusion  of  the  signal  in 
the  matrix  inversion  and  subsequent  mismatch  has  been  treated  as  a  special  case  of  inter¬ 
ference  rejection.  A  simple  expression  for  this  anomalous  signal  suppression  has  been 
presented.  This  signal  suppression  can  only  be  significant  if  the  signal-to -noise  ratio 
would  be  large  at  the  output  of  an  optimum  beamformer  steered  perfectly  in  the  signal 
direction. 

Various  beamformer  optimization  problems  have  been  considered.  Constraining  p 
was  found  to  be  equivalent  to  adding  a  spatially  incoherent  component  to  the  signal  field 
prior  to  unconstrained  optimization.  Constraining  was  found  to  be  equivalent  to  adding 
a  spatially  white  component  to  the  noise  field. 

In  many  important  situations  the  simple  beamformer  structure  is  not  optimum.  A 
more  general  array  processing  configuration  has  been  presented  which  provides  for 
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incoherent  combination  of  the  outputs  of  a  number  of  simple  beamformers.  Optimization 
problems  associated  with  this  processor  have  been  solved  providing  generalization  of  the 
results  for  simple  beamformers. 
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Figure  1.  General  Beamformer  Configuration. 
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Figure  2.  General  Array  Processor. 


